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We investigate entanglement transfer from a system of two spin-entangled electron-hole pairs, 
each placed in a separate single mode cavity, to the photons emitted during their recombination 
process. Dipole selection rules and a splitting between the light-hole and the heavy-hole subbands 
are the crucial ingredients establishing a one-to-one correspondence between electron spins and 
circular photon polarizations. To account for the measurement of the photons as well as dephasing 
effects, we choose a stochastic Schrodinger equation and a conditional master equation approach, 
respectively. The influence of interactions with the environment as well as asymmetries in the 
coherent couplings on the photon-entanglement is analyzed for two concrete measurement schemes. 
The first one is designed to violate the Clauser-Horne-Shimony-Holt (CHSH) inequality, while the 
second one employs the visibility of interference fringes to prove the entanglement of the photons. 
Because of the spatial separation of the entangled electronic system over two quantum dots, a 
successful verification of entangled photons emitted by this system would imply the detection of 
nonlocal spin-entanglement of massive particles in a solid state structure. 
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I. INTRODUCTION 

The macroscopic scalability of its architecture is 
among the fundamental requirements which are to 
be met by any sound implementation of a quantum 
computer^. The ability of transferring entanglement 
from spatially stationary building blocks, e.g. condensed 
matter embedded storage registers, to flying qubits, e.g. 
photons, would be an important step towards building 
such a computer. In this scenario, quantum communica- 
tion via photons would replace the concept of data buses 
approved in the framework of classical computing. The 
main task of our work is to extensively investigate such 
a transfer process including the theoretical description 
of different kinds of non-idealities and measurements. 

In the past, several proposals concerning the pro- 
duction of entangled photons employed the indis- 
tinguishabilit y of two d ecay paths within a biex- 
citon cascad<P3EEEEEED During the last years, 

gre at experimental pr ogress has been made in this 
fiel( jiu | ii | i2 | i3 | i4 | i5 | ib | i7 | In these systemS! entanglement 

is generated during a coherent twofold decay process, 
whereas we intend to transfer entanglement in a con- 
trolled way from electron spins to photon polarizations. 
The functionality of a similar device performing the lat- 
ter task has been proposed and studied in Ref. [TS] using 
a master equation approach. However, a shortcoming of 
this proposal has been the necessity of postprocessing 
steps to disentangle the electronic system from the pho- 
tons. Our system avoids this additional effort by using 
an entangled hole pair and an entangled electron pair 
instead of assuming only the electrons to be entangled 
(for similar proposals see Refs. I19l20p . In this work, 



we go beyond existing literature and investigate the 
influence of dephasing in the electronic system on the 
entanglement of the emitted photons with the help of 
a quantum trajectory pictured This formally involved 
treatment in the framework of open quantum systems^! 
is justified by showing unambiguously that some in- 
teresting features of the system's behaviour as to the 
production of entangled photons cannot be understood 
properly without an unravelling of the quantum master 
equation in terms of quantum trajectories. 

In the following, we present a schematic of our entan- 
gler's functionalit y. A double dot device consists of a lat- 
eral quantum dot ! * 24 * provided with an electrostatically 
tunable constriction potential used to divide the single 
dot into two separate quantum dots^. We assume one 
double dot to be charged with two electrons and another 
one with two heavy holes (HH) . Both pairs of charge car- 
riers are then supposed to relax into their groundstate 
which is a spin singlet state. By turning up the constric- 
tion potential in each of the double dots we separate the 
initially indistinguishable particles orbitally to yield two 
entangled spin pairs. Thereafter, respectively one elec- 
tron and one hole are transported to an optical quantum 
dot via nanowireiPfil so that the spin pairs finally furnish 
two entangled exciton states delocalized over two spa- 
tially separated optical dots A and B each surrounded 
by a single mode cavity (see Fig. [T]). It will be shown 
that by optical recombination of the two entangled ex- 
citons a polarization entangled two photon state of the 
type 

I*/) = |0) e - h+ ® ^(k>B> + (i) 
can be produced which is a Bell state and thus maximally 
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Figure 1: Schematic of the photon entangler. The neigh- 
bouring circles denote lateral single dots gained by dividing 
double dots. The optical dots (pink boxes) contained in the 
nanowires (white stripes) are distinguished by naming them 
A and B. The black dots depict electrons, the small circles 
HH. Photon cavities are shown in blue oval areas around the 
optical dots. 



entangled. The two circularly polarized photons are 
then supposed to be measured after they have leaked out 
of the cavity surrounding dot A and dot B respectively. 
Preliminary experimental progress along these lines 
has been reported in the context of spin light-emitting 
diodes™ 

This article is organized as follows: In Section |TTJ a 
consistent modelling of the quantum dot photon entan- 
gler and its dynamics in the framework of open quantum 
systems will be presented. Within this model the photon 
entangler is assumed to be coupled weakly to its environ- 
ment which exhibits the basic properties of a heat bath so 
that the Born-Markov approximation 29 is applicable. A 
quantum stochastic differential equation (QSDE) will be 
derived as the equation of motion of our system of inter- 
est. The solutions of this equation are the quantum tra- 
jectories representing the time evolution of the reduced 
system's pure state. Furthermore, we introduce a level 
of description which is realistic as measured by the infor- 
mation a conscious observer measuring the emitted pho- 
tons gains about the system state. The relevant e quation 
of motion is then a conditional master equation 1 4 * 21 * 30 ^. 
In Section |III[ we discuss our simulation results for the 
quantification of entanglement which becomes manifest 
in the violation of the CHSH inequalitjEH and the visibil- 
ity of interference fringes proving the existence of quan- 
tum coherence in a two photon Hilbert space. We work 
out why a treatment beyond the quantum master equa- 
tion approach is adequate by capturing the fingerprints 
of the quantum trajectory picture in the conditional den- 
sity matrix the conscious observer is aware of. Section |TV| 
is aimed to sum up the most important findings and to 



give an outlook as to possible future work built on this 
article. 



II. MODEL & METHODS 

A. Initial state preparation 

To be capable of transferring entanglement from the 
electron hole excitations in the optical dots to the photon 
pair emitted during their recombination, the optical dots 
must be charged with an appropriate initial state. This 
may well be accomplished by taking the following steps: 

• Charge each lateral quantum dot with a singlet of 
two elect rons/HH. 

• Divide each lateral dot coherently into two (as 
shown in Fig. [I]) by turning up a constriction po- 
tential. 

• Get the charge carriers transported into the optical 
quantum dots. 

• End up with two entangled exciton states delocal- 
ized over the two optical dots A and B. 

If all these steps are taken in an absolutely coherent way, 
the spin state of the two e"-HH excitations in the optical 
dots will yield 

1*0 =l(\n)A\m B -\n)A\ii\)B - \utu\n)B+ 



dark 
dark 



(2) 



where the double-arrows denote the HH-hole spins and 
the simple ones represent the electrons (see Fig. [2|. 



B. Entanglement transfer due to optical dipole 
transitions 

Given the initial state of Eq. ([2 1 our device is now 
ready to perform its key-task, namely the entanglement 
transfer from the electronic system to a photon pair. Due 
to the dipole selection rules A J — ±1, Amj = 0,±1 
only certain recombinations of e~-HH pairs are allowed 
(see Fig. [21). The one particle states of the e~-HH ex- 
citations which are not allowed to recombine differ by 
2 units of H in their spin z-projection. These not radia- 
tively active states are called dark exciton states (see Eq. 
(J2|). The cavity mode for both subsystems, A and B, is 
assumed to be on resonance with the electron-HH re- 
combination freqency Uq and to propagate in z-direction. 
The z-direction is defined by the quantization axis of an- 
gular momentum which arises due to HH-light hole (LH) 
splitting. This typical splitting comes from the confine- 
ment of charge carriers in the optical dots and is am- 
plified by mechanical strain to which the quantum dots 
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|T) m s = + \ I) rn s = -\ 




Figure 2: Dipole allowed recombinations establishing a one- 
to-one correspondence between photon helicity and electron 
spin direction, mj = +| denotes a HH-state in which a +|- 
electron is missing. The light hole (LH) levels are assumed to 
be far off-resonant as measured by the cavity linewidth. 



are often exposed during the growth process. For pho- 
tons emitted in z-direction, the only direction allowed by 
the cavity, the polarization will be circular which dis- 
plays the conservation of angular momentum along the 
z-axis and implies Amj — ±1. In other directions the 
radiation of elliptically polarized light would be possible 
which precludes the desired one-to-one correspondence 
between photon helicity and electron spin projection^. 
The necessity of the HH-LH splitting becomes clear look- 
ing at Fig. [2] because the correspondence between circu- 
lar photon polarization and electron spin would also be 
destroyed by transitions involving LH states. These tran- 
sitions are prevented by the assumption that our cavity 
provides only a single radiation mode on resonance with 
the e"-HH energy gap. However, it should be mentioned 
that the dark exciton states limit the fundamental effi- 
ciency of our entangler to 50%. 

Accounting for the selection rules and omitting the dark 
states our device will produce a maximally entangled two 
photon state (see Eq. [T]), as its final state after recombi- 
nation in both optical dots A and B. The dots are left in 
the vacuum state, i.e. without any electronic excitations 
and are therefore not entangled with the photons. 



1. Hilbert space for the two-dot system 



We model the conduction band and the HH band with 
a single level for each spin projection assuming that 
other subbands are shifted far away, as measured by the 
linewidth of our cavity mode, due to spatial confinement. 
Furthermore, we assume that the energy levels are iden- 
tical in both subsystems A and B. An overview over all 
states reachable from the initial state Eq. (pi) is given by 
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States with one excitation 
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Vacuum state without any excitation 

|Q) = \0)a ® |0)b ® |0), h . (3) 



C. System states and dynamics 

This section is intended to state more precisely the 
behaviour of our system. The system's Hilbert space and 
its coherent dynamics as well as models for environmental 
influences are introduced. For clarity, we start with a 
subsystem containing only the two optical dots before its 
coupling to the lateral dots provided by the nanowires is 
investigated. 



The blocks are arranged by number of excitations. 
a^ k , i = ±, k = A, B creates a photon with polarization 

i and frequency loq in cavity k. h»k , k = A, B creates 
a heavy hole with spin down and annihilates an electron 
with m.j = + 1 in the valence band of dot k respectively. 
c| fc , k = A,B creates an electron with mj = +^ in the 
conduction band of dot k. Double arrows denote HH, sin- 
gle arrows denote electrons. Note that the dark states, 
the dynamics of which is irrelevant as to the production 
of photons, are neglected. 
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2. Coherent dynamics 

The Hamiltonian of the two-dot system seen as a closed 
system is defined by the interaction of the electronic de- 
grees of freedom with the radiation field which is, given 
a single mode cavity on resonance surrounding each op- 
tical dot, governed by the Jaynes-Cummings modeP^I. In 
the interaction picture, i.e. after subtracting the nonin- 
teracting Hamiltonian which is nothing but the sum of 
the excitation number operators, the Hamiltonian for one 
optical dot reads in the rotating wave approximation 

iff = q k c^ k h^ k (g>a a k+v k c^h^f&a^+h.c; k = A,B; 

(4) 

with the cavity coupling strengths {qi,Vi} and the an- 
nihilation operator of the cavity mode a CT i_ of cavity i 
with polarization er + and cr_ respectively. The total 
interaction Hamiltonian of both dots is then given by 

H = H[ A) ® l( B ) + l^® H[ B) . (5) 

The interaction Hamiltonian commutes with the total ex- 
citation number. The coherent unfolding is therefore not 
capable of connecting states of different excitation num- 
bers (see (Eq. [3])). This kind of transitions can only 
be fulfilled by incoherent processes which will be treated 
next. Note furthermore that the coherent dynamics gen- 
erated by H factorizes so that the total time evolution 
operator can be written as a tensor product, that is 

U(t,to) = Uf(tM)®U?{t,t a ) = e - ,; < ( *- to) (g)e-^ B(t - to) . 

The following analysis will therefore be done in the sub- 
space of one optical dot. 



3. Dissipation channels 

Now we account for the imperfection of the cavity, 
which is essential for later purposes because we are only 
capable of measuring photons which have leaked out of 
the cavity. The cavity mode couples to the surround- 
ing radiation field which is assumed to be initially in the 
vacuum state, i.e. in thermodynamic equilibrium at zero 
temperature. The corresponding dynamics is specified 
by the following ansatz for the interaction picture leak- 
out Hamiltonian: 

H£(t) = [a\ h+ ® B+(t) + +, ® B^(t)}+h.c, (7) 

where k = A,B and ujq denotes the resonance frequency 
of the electron-HH recombination which is equal to the 
frequency of the cavity mode. The bath operators B~£ 
are given by 

Bt(t) = (J2^- iUjtb 3 + ) (8) 
i 



with the leakout coupling strength £j of the cavity mode 
to mode j of the external field. For a definition of B^ one 
needs to exchange + and - in the above definition of Bf . 

The 6 . ' are the annihilation operators of the radiation 
field mode j outside the cavity k = A,B with polarization 
a±. We approximate the radiation field modes {bj} out- 
side the cavities with a continuous spectrum with DOS 
D(lo). Performing the Markov approximation we can cal- 
culate the jump operators {Ai} and their respective rates 
{Pi} relevant for the reduced system's dyn amics in the 
framework of continuous measurement theorjl21122J Up to 
second order perturbation theory in Hj := + Hi if), 
we get 
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(9) 



7 = 27rD(a;o)|^(wo)| 2 is the resonant cavity leakout rate. 
The dynamics of the two-dot system including the cavity 
leakout dissipation channels is depicted in Fig. [3] 




|10)«^|9) |12)<^|11> |14)A|13) |16)*^+|15) 




Figure 3: Open quantum system dynamics of the two-dot sys- 
tem. The rows are arranged by excitation number. The blue, 
down pointing arrows represent photon measurements. Co- 
herent couplings are drawn with black arrows. The dephasing 
operators are not pictured for clarity. 



Up to this point the considered dynamics is rather ide- 
alized because it does not account for any uncontrollable 
interaction which is able to dephase the delocalized sin- 
glet states into product states. Such an interaction will 
now be introduced by adding four more dissipation chan- 
nels which distinguish the spin polarizations of the dif- 
ferent electron- HH excitations. This model is equivalent 
to an independent dephasing of the four single particle 
spins since the spins occur only as part of an e~-HH ex- 
citation in our dynamics. The projectors on the relevant 



5 



exciton numbers, i.e. our dephasing operators read: 
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:= h^ B h^BC^ B c^B 


(10) 



Note that if coherences between the two excitons oc- 
curring in the initial state get lost, neither the elec- 
trons nor the holes will be entangled any longer. This 
would also preclude the entanglement of the emitted 
photons, of course. The dephasing jump operators 
™udi ^dui^udi ™du are assumed to act at a rate fli = k, i — 
5, . . . , 8 which we call the dephasing strength. These op- 
erators also act only on the subsystem Hilbert spaces 
T~(-a,T~(-b of the two-dot system. 



4- Stochastic equations of motion 

Accounting for all dissipation channels just introduced, 
the quantum stochastic differential equation (QSDE) for 
the two-dot system is given by 



d*(t) = -learnt + J2 ( pffjfj - *w) dNk (*)' 

with the nonlinear deterministic generator 



k=l 



k=l 



= :H 



(12) 

The dissipation channels governed by the dephasing jump 
operators are of course not controllable for a conscious 
observer so that, if consideration is given to a realistic de- 
scription, the pure state dynamics specified by the latter 
QSDE is no longer adequate in the presence of dephasing. 
A level of description, containing exactly the knowledge 
a conscious observer may gain by perfectly detecting the 
leakout photons, is provided by the following conditional 
master equation^, 



dp s = -i[H,p s ]dt 

4 

E 



A kPs A\ _\ dNk + lTi . {A ^ Akps}psdt 



fc=i k \Tr{AtA kPs } , 

{A{A k , Ps }dt} + « ^A kPs A\ - l -{A\A klPs }^ dt. 

(13) 



7 i 
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Setting the stochastic increments for the photon detec- 
tions {dNk] to zero in Eq. (13 1, we get 



dp. 



, t = -i[H, Pc ] + lTv{A{A kPc } Pc 



k=l 
8 



fc=5 



k=l 



(14) 



This evolution, conditional on the system not emitting 
a photon out of the cavity, takes place between the two 
photon detections and before the first one. It has been 
showrPl that the nonlinear term in the latter equation 
is only relevant for the normalization of the conditional 
density matrix. If we renormalize the density matrix at 
any given time by hand before calculating the expecta- 
tion value of any observable it is thus sufficient to solve 
the following linear ordinary differential equation: 



-i[H, Pc ] + J2 *A k pcA{ - £ f {A{ A k , Pc }, 



dt 



fc=5 



k = l 



where Pc is the unnormalized conditional density matrix. 



5. Coupling between optical and lateral dots 

Concerning the transport of charge carriers from the 
lateral to the optical quantum dots our system can again 
be decoupled into two parts called A-side and B-side (see 
Fig. [TJ.The transport of electrons and HH through the 
nanowires is modelled by a simple tunneling process. Dif- 
ferent tunneling rates for the different spin projections 
can be chosen to account for asymmetries of the meso- 
scopic wires. An A-B-dependence of the tunnel couplings 
models imperfections in the production of identical com- 
ponents. The lateral dots are treated as one level sys- 
tems for each spin state. A diagram of the relevant states 
when consideration is given to the coupling of lateral and 
optical dots can be found in Fig. [3] This description 




Figure 4: Relevant states for the transport of electrons 
and HH between optical and lateral quantum dots, k = 
A,B. {s'y,,} denote the tunneling rates for up/down elec- 
trons, {t^/d} are the respective couplings for the HH. 
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introduces eight new parameters, the four tunnel rates 
{s^i = u,d; k = A, B} for electrons and the four respec- 
tive couplings {t*\i — u,d; k = A, B} for HH. 
The dynamics of our system is assumed to start hav- 
ing the particles localized in the lateral dots while the 
spin-state itself is not altered compared to Eq. |2]). The 
Hamiltonian H$ for the coherent dynamics of the whole 
system reads 



Ha- H 4-{ \" c fc J r (L) 4- s fc r f r [L) 4- t k h {L) 



k=A,B 



td, Mj,fc h M), + ft-.C.}, 



(16) 



where the index (L) denotes the lateral quantum dots. 
In the simulations for the coupled system we neglect de- 
phasing so that we don't need to expand our dcphasing 
model to the enlarged system. The difference between 
the whole system and the two-dot system is then simply 
the substitution of H s for H in the respective stochas- 
tic equations of motion since the cavity leakout concerns 
only the optical dots anyway. 



D. CHSH violation 

Now, we prove and quantify the capability of our sys- 
tem to produce entangled photons by simulating a CHSH 
inequality 31 violating measurement. We focus on basic 
physical features which are essential to understand the 
advantages of the quantum trajectory picture. For sim- 
plicity we assume that the two photons are measured 
simultaneously, which could experimentally be realized 
by a post selection rule only keeping track of coinciden- 
tal measurements. Furthermore we only investigate the 
two-dot system not describing the coupling to the lateral 
dots. 

Without dephasing, our system is supposed to produce 
the maximally entangled two photon state \^ p h) — 



Wa°b)) ( see E q- 



^\\ (t a (X b) + \°a g b)) l see fc q- An appropri- 

ate measurement scheme to verify the entanglement of 
the state measured in the simulation should accordingly 
be designed such that the CHSH inequality will be vio- 
lated maximally by \^ p h), that is B(^> p h) = (QS)v ph + 
(RS)* ph + (RT)^ ph - (QT)^ ph = 2V2. This condition is 
met by the subsequent set of observables. 



Q = 

s = 



Za, 

—Zb — Xb 



T = 



R=-X A 
Zb — Xb 



(17) 



Note that Q and R are observables of subsystem A while 
S and T are observables of subsystem B. The combina- 
tions QS, QT, RS and RT are each measured with prob- 
ability 0.25 in all simulations. 

To calculate the expectation values of the relevant ob- 
servables we have to solve the stochastic equation of mo- 
tion, Eq. (13) and Eq. (11) respectively to obtain the 



We employ several methods to accomplish this integra- 
tion: 

A first choice providing an exact solution is the in- 
tegration of Eq. ( 15 1 to calculate the density matrix 

Mil 



Pc(t) 



Ti'Pc(r) 

ues are then given by 



The relevant desired expectation val- 



(XY) 



Pc(T) 



Tr {p c (r)XY} 
Tr{p c (r)} : 



X = Q 7 R; Y = S,T. 

(18) 

A second method to obtain these quantities is the 
numerical integration of the QSDE, Eq. (Ill, by creating 



a stochastic ensemble of quantum trajectories (see Refs. 
1211221 for a detailed discussion). The expectation values 
of the mentioned observables are obtained determining 
the respective measurement results for each trajectory 
and then calculating the ensemble average. This ap- 
proach reproduces the exact solution when the ensemble 
size goes to infinity. 

A third technique is to calculate the integral represen- 
tation of the propagator of the underlying piecewise 
deterministic process PDP with the help of a Hilbert 
space path integral. By this we would calculate 
the statistical contribution of any plausible quantum 
trajectory and therefore obtain the same results as 
when averaging over an ensemble of infinite size. To 
realize the computation of the path integral we have to 
introduce a cutoff c which leads to the neglect of all 
trajectories containing more than c dephasing jumps. 
This approximation is reasonable since we are in a weak 
coupling regime where the Born-Markov approximation 
is applicable. Hence, we can think of k as small parame- 
ter which suppresses trajectories with i dephasing jumps 
with a prefactor of k 1 . The path integral technique then 
converges against the exact solution as c goes to infinity. 



E. Quantification of entanglement with a single 
interference setup 

Our second approach makes use of the fact that the 
dimension of the polarization Hilbert space tO^g of 
the two leakout photons is reduced in virtue of opti- 
cal selection rules. They forbid the emission of photon 
pairs with equal circular polarization which cuts down 
U 



ph 
AB 



to a two dimensional Hilbertspace Ti^g spanned 



reduced system's state at the photon measuring time r. 



by the basis M = {Ic^cr^), W^crg)}. The task of 
investigating polarization entanglement is then simply 
to distinguish a coherent superposition like a\a^ag) + 
f3\(j\cTg), I a 1 2 + \(3 1 2 = 1 from an incoherent mixture 

|a|Vi CT BXfA°sl + l/ ? l 2 lf r A cr B>( cr A cr sl and t0 deter- 
mine the ratio of the amplitudes a and (3. These proper- 
ties can be characterized with the help of a single inter- 
ference experiment. 

Our interference setup is closely related to that proposed 
in Ref. |U For later convenience we prepend a device 
which flips the polarization of photon B (see Fig.|5|, that 
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PBS PBS 




PBS 



Figure 5: Interference setup for the postprocessing of photons 
emitted by subsystem k = A and k = B respectively. PBS 
stands for polarization beam splitter, the half wave plate is 
denoted with HWP. Th performs a polarization flip on sub- 
system B leaving subsystem A unchanged. 



looking at the following equation: 




The one-to-one correspondence between entanglement 
and visibility becomes manifest in the subsequent rela- 
tion between V and the Bell parameter B max of a mea- 
surement scheme maximizing B for a given input stated 



B max = V2(l + V). (25) 
III. RESULTS & DISCUSSION 



is a unitary operation 



A. Decay of entanglement due to dephasing 



T = T A ® Fb = 1 



ph 



(19) 



to the interferometer which without T performs the 
transformation 
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(20) 



The effect of the whole setup shown in Fig. [5] on the vec- 
tor of photon annihilation operators a :— (a a A : a a A) T £3 

(a a B 1 a cy B) T reads 



a' = OoT(a). 



(21) 



Note that the unitary operation O o T can be factor- 
ized to a tensor product of two local operations in the 
respective one photon Hilbert spaces. This is of great 
importance, because otherwise it could alter the entan- 
glement properties of the input state. 
Following the analysis in Ref. 4, we define the visibility 
of interference fringes V by 



V(t B ,t A ) = 



2\Z(t B ,t A )\ 



X{t B ,t A )+Y{t B ,t A )' 
where the quantities X, Y, Z read in our case 
X = (: h a B(t B )h a A(t A ) :) = 

= Tr{<Vf2"is,U \ a a AT tA ,o{p a {Q)) a\^j a^} 
Y = (: h a B(t B )h a A(t A ) :) = 

= Tr {a a B% B<tA {a a AT tAfl {p s {Q))al^j a^j 



(22) 



(23) 



T. . is the conditional time evolution operator represent- 



ing the flow of Eq. (14 1, n a h, i = ±, k = A,B are the 



photon number operators of the respective cavity modes 
and the colons denote normal ordering. The name visi- 
bility of interference fringes for V is immediately justified 



In order to understand the basic phenomenology of 
entanglement decay due to dephasing we set Uj = qi = 
1, i = A, B and j = 0.1 so that the coherent coupling 
of the electron-HH excitations to the cavity mode is by 
one order of magnitude stronger than the leakout rate 
and we have total symmetry between subsystems A and 
B. The influence of asymmetries in the model parameters 
is investigated with the help of the interference approach 
in Section IIIB All energies arc measured in units of 



q B throughout this work. The dimcnsionlcss timcscalc is 
then fixed by the inverse of this energy- unit. 



1. Results for the CHSH violation 

The dephasing strength k is varied ranging from k = 
to k = 0.2. Without dephasing (k — 0) the expected 
value of B = 2\/2 can be reproduced for any given emis- 
sion time. For a totally dephased product state, B = y2 
holds for our measurement scheme presented in Section 
|IID| Hence, we expect in the presence of dephasing an 

asymptotic behaviour like B(t) y/2. 
The time dependence B(t) of the Bell parameter for 
k = 0.05 calculated with the three different integration 
techniques described in Section |IID| is shown in Fig. [6] 
The concurrence of the three plots is quite well, showing 
no systematic deviations between the numerically exact 
solution (red solid) and the noisy stochastic ensemble av- 
erage (blue dotted). The path integral calculation with 
cutoff c = 2 (green dashed) is almost congruent with the 
numerically exact solution for small emission times and 
shows significant aberrations systematically overestimat- 
ing the entanglement only for t > 10. That is because 
the weight of trajectories which are of higher order in 
k becomes significant for greater emission times so that 
we would need a higher cutoff to obtain adequate ac- 
curacy. All plots pictured in Fig. [6] have two remark- 
able features. First, a net decay of B due to dephas- 
ing, enveloping the oscillations can be found. Second, 
oscillations with the same frequency as those of the co- 
herent emission probability (see Fig. [T]) which is pro- 
portional to the population of the two photon states, i.e. 
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plotted for k = 0.05 calculated by solution 
with a Hilbert space path integral with cutoff 
by solution of Eq. ( 15 1 (red solid) and 



Figure 6: B(t) 
of Eq. {IT} 

c — 2 (green dashed) 
by averaging over an ensemble of 2000 quantum trajectories 
(blue dotted). 



V{t) ~ (|7)(7| + |8)(8|) \%(t)) , are visible. The un- 




Figure 8: Time dependence of the two photon population 
P for a totally coherent trajectory (blue dashed) and a 
trajectory with one dephasing jump (red solid). The red 
curve shows clearly, that the phase coherence of the Rabi- 
oscillations in the blue plot gets lost due to dephasing. 




Figure 7: Time dependence of the two photon emission prob- 
ability V(t) for k — 0.05. The distribution is not normal- 
ized, since only simultaneously emitted photons are consid- 
ered. The plot displays Rabi-oscillations generated by the 
Jaynes-Cummings type coherent dynamics. 

derstanding of this complex time dependence lies at the 
heart of our quantum trajectory description and is one 
of the most interesting results of this work. The Jaynes- 
Cummings type dynamics brings about Rabi-oscillations 
between exciton and photon states in each of the optical 
dots which penetrate to the photon emission probabil- 
ity V (see Fig. j7|). At the minima of V the coher- 
ent dynamics forbids the emission of a photon because 
only the electron-HH states are populated and no pho- 
tons are present. A dephasing event before the supposed 
measurement randomizes the phase of the coherent Rabi- 
oscillations. This enables trajectories which have been 
exposed to dephasing to emit photons at any instance 
in time, in particular at times where the coherent Rabi- 
oscillations enforced a purely excitonic state (compare 
the two plots in Fig. [8}. For these measuring times 
t rn the dephased trajectories contribute to a very large 



amount, up to hundred percent at the sharp minima of 
V, which explains the minima of B(t) with B(t m ) w V2 
at these points. 

The net decay of B(t) is quite intuitive since the probabil- 
ity that a dephasing jump occurs increases monotonously 
with time. To investigate the long-time behaviour of the 
entanglement decay it is obviously (see Fig. [9} accurate 
to fit an exponential function to the net decay of en- 
tanglement. The dependence of the exponential decay's 




Figure 9: Long-time decay of B(t) plotted for « = 
0.075 (red solid) and fit) = V2(l + exp(-0.585i)) (blue 
dashed). The blue exponentially decreasing curve envelops 
the net decay of entanglement. 

measure on the dephasing strength is studied by plotting 
the time t v at which B(t v ) — 2 holds, that is to say at 
which the CHSH inequality is met, versus the inverse de- 
phasing strength (see Fig. 10 1. t„ increases linearly with 
1/k which indicates that the measure of the exponential 
net decay is reciprocally proportional to the dephasing 
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50 



1/k 



Figure 10: The time t v at which the CHSH violation vanishes 
plotted against - (blue points). The concurrence with the 
linear fit t(l/n) = -0.3166 + 1.1529^ (red line) is on the spot. 




strength n. For 1/k — ► 0, t v goes up to a small error 
to zero, which displays infinitely fast entanglement decay 
for k — ► oo as one would intuitively expect. 



2. Results for the interference approach 

The decay of entanglement as measured by the visibil- 
ity of interference fringes V(t B ,t A ) is shown in Figure 
[TT]for different dephasing strengths n = 0.05, 0.10, 0.20, 
7 = 0.1 and Vi — q .j = 1.0, calculated by numerical inte- 



gration of Eq. (15). The behaviour is a straightforward 



generalization of the phenomenology discussed in Section 
|III A 1| to nonsimultaneous photon measurements. We 
again see a net decay with increasing emission times and 
characteristic oscillations following the Rabi-oscillations 
of the emission probability (see Fig. 12 1. The emission 
probability V is defined in a straightforward way by 

V(t B ,t A ) = 

Tr {^f Tt B ,t A [KfPitA^A | h a B I Tr{p(t A )} 



7 2 (*+n 



.1.0. g. t A < t B , 



(26) 



where X,Y are the quantities defined in Eq. (23 1 but 
with the unnormalized density matrix p introduced in 
Eq. ( 15 1 used to calculate expectation values. Tr{p(iyi)} 



is the probability that no photon emission occurs before 
t A - Tt B ,t A denotes the time evolution operator governed 



by the flow of Eq. ( 15 1 and brings about the conditional 
time evolution of the nonnormalized density matrix p 
between t A and t B - As analyzed in the previous sec- 
tion we assume the oscillations of the visibility to disturb 
the functionality of the photon entangler only minimally, 
since the emission probability in the valleys of V(t A , t B ) 
is minimal. The major contribution to the mean value of 
the visibility, 



V = 



dt A / dt B V(t B ,t A )V(t B ,t A ), 



(27) 





Figure 11: Decay of V(tB,tA) for k = 0.05 (top), k = 
0.10 (middle) and k = 0.20 (bottom), calculated by solution 
of Eq. 1 15 1. 7 = 0.1 and Vi — qi = 1.0 in all plots. 



is thus expected to be given by the net decay of V. 
To confirm this hypothesis, we calculate and plot the 
visibility for a small dephasing strength k = 0.01 for 
which the net decay of entanglement is very slow (see 
Fig. 13 1. However, the oscillations bringing about 



distinctive valleys in the visibility suggest that the mean 
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Figure 12: V{t B ,t A ) for 7 = 0.1 (top), K - 0.05 and v t 
qi — 1.0, calculated by solution of Eq. (15 1. 



1. Asymmteries in the two-dot system 

We distinguish two types of imbalances, the first of 
which is an imbalance between A and B whereas the 
second asymmetry is due to unequal couplings for the 
different polarizations. A pure A-B-asymmetry, that is 
to say VA = q_A^VB = qB, does not affect the en- 
tanglement of the emitted photons at all. That is be- 
cause it does neither break the symmetry between the 
two kets Ic^Og), \ a A a ^)i t ne coherence between which 
makes up the entanglement, nor introduce any asym- 
metry between the two polarizations on single photon 
space. If we set up a polarization asymmetry by setting 
Va = 1-1, qA = 1.0, vb = Qb = 1-0 only on subsys- 
tem A, this symmetry is broken causing a complex in- 
dependence of V whereas the visibility is independent of 
ts (see Figure 14). The emission probability V shows 




Figure 13: Decay of V(tB,t A ) for k, = 0.01, calculated by 
solution of Eq. 7 = 0.1 and Vi = qi = 1.0 . The net 

decay is very slow but distinctive valleys are still visible for 
this weak dephasing strength. 



entanglement of the emitted photons could be signifi- 
cantly reduced. The mean visibility for this parameter 
set yields V ~ 92.1% which justifies the statement that 
the oscillations do not disturb the functionality of the 
entangler to a large amount, because the loss of about 
eight percent is almost entirely covered by the weak net 
decay. 



B. Influence of asymmetries in the coherent 
couplings 

Now, we focus on the discussion of a different kind 
of non-ideality, namely the influence of asymmetries in 
the Jaynes-Cummings couplings {qi, Vi\i — A, B} on the 
visibility V . From now on we neglect dephasing, i.e. we 
set k = 0. 




Figure 14: V{t B ,t A ) (top) and V{t B ,t A ) (bottom) for v A = 
1.1, q A — 1.0, vb — qs — 1-0 , calculated by solution of Eq. 
( 15 1. 7 = 0.1 and k = in both plots. 



the regular Rabi-oscillations dependent on ts whereas 
its i^-dependence is more irregular, governed by the com- 
peting oscillation frequencies va, qA (see Figure 14). 

If a polarization asymmetry is introduced in the same 
way on both subsystems, for example by setting va — 
vb = 1.1, qA = qB = 1.0, the symmetry in the two exci- 
tation subspace of our Hilbert space (i.e. the symmetry 
between left and right hand side of the first row in Fig. 



11 




by considering the coupled Hamiltonian H$ (see Eq. 
(16l) as the generator of our coherent dynamics. In- 



Figure 15: V{tB,tA) (top) and V{tB,tA) (bottom) for va = 
vb = 1.1, q_A — q_B — 1.0, calculated by solution of Eq. (15 1 
7 = 0.1 and k = in both plots. 



[3]) is not broken. However, the obvious asymmetry in 
the one exciton subspace will cause an effect as soon as 
one photon is emitted. If the two photons are emitted 
simultaneously V = 1 thus holds (see the diagonal of the 



upper plot in Fig. 15 1, while the behaviour of the visibil- 



ity becomes rather irregular for arbitrary emission times 
(tA, tg), nevertheless preserving t^-is-symmetry (see the 
upper plot in Figure 151. The emission probability (see 
the lower plot in Fig. 15 1 exhibits oscillations with com- 
peting frequencies Vi = 1.1, qi = 1.0 in both directions 
tA , is also keeping the A-B-symmctry. To sum up the 
latter analysis as an advice to the experimentalists, we 
state that the cr+-cr_-symmetry of the cavity couplings 
should be the preliminary goal as to maximize the vis- 
ibility. Furthermore, if the requirement of polarization- 
independence can not be met but A-B-symmetry is pre- 
served, then the post selection rule t& — tg can guaran- 
tee perfectly entangled photons within the framework of 
our model. 



2. Asymmetries in the tunnel couplings 

In this part of the discussion we extend the coher- 
ent dynamics to the whole system's Hilbert space. That 
means we include the tunnel coupling to the lateral dots 



vestigating the influence of asymmetries in the tunnel 
coupling strengths {sf,tf|i = u, d; j = A, B} we ob- 
serve a phenomenology similar to that occurring when 
consideration is given to imbalanced cavity couplings 
{i)i,qi\i = A, B} within the two-dot system. The idea 
of correcting the effect of these imbalances by system- 
atically generating asymmetries in the tunnel couplings 
is therefore straightforward. Unfortunately, such a cor- 
rection as measured by an improvement in mean visibil- 
ity could not be achieved varying the relevant model pa- 
rameters. Given symmetrical parameters for the two-dot 
system, the visibility again is inert regarding asymme- 
tries which break A-B-symmetry but preserve electron- 
HH and polarization symmetry, that is s A = s d — t„ = 



t A 



if 



t B 



The visibility also stays 
at its maximal value V = 1 if we break the electron- 
HH symmetry by setting s k = s k d ^ t* = t d , k = A, B. 
This observation displays the fact that we have to intro- 
duce an asymmetry which is polarization selective, if we 
want to touch the polarization entanglement of the emit- 
ted photons. The simplest way to do that is by choosing 
s u s d — = on one dot k- The resulting visibility 
V(tB,tA) and emission probability 7 5 (t_B, tyi) are shown 
in Figure |16| As expected, the visibility is only depen- 
dent on tA , since the polarization asymmetry concerns 
only subsystem A. The emission probability exhibits the 
same behaviour as in the symmetrical case concerning its 
i^-dependence, whereas its t^-dependence is influenced 
by the two competing coupling frequencies s A ^ s^. 
Let us now create a polarization asymmetry which pre- 
serves A-B-symmetry by setting s£ 7^ s d — t\ = t k d 
on both dots k = A,B. Like for the two-dot system we 
then expect the coherence of our relevant states to be 
untouched as long as both excitations are present in the 
system. The asymmetry becomes only relevant after the 
first photon has been measured. The visibility should 
therefore be maximal for tA — tg- This behaviour, to- 
gether with the expected A-B-symmetry of V, V can be 
seen in Figure [17] 
The concluding suggestion to the experimentalist striv- 
ing for the measurement of entangled photons is very 
similar to the discussion of asymmetries in the Jaynes- 
Cummings couplings but now with the tunnel couplings 
{s k ,t k } in the roll of {v kl q k }. 



IV. CONCLUSIONS 

This final section is aimed to review the most impor- 
tant findings and to give an outlook as to future research 
which may be built on this work. 

The first remarkable observation is the appearance of 
characteristic oscillations in the emission time dependent 
decay of two photon entanglement which are congruent 
with the Rabi-oscillations of the emission probability 
V . The understanding of this phenomenon has required 
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Figure 16: V(t B ,t A ) (top) and V(t B ,t A ) (bottom) for s£ = 
1.1 ^ sf = t£ = tj: = sf = tf = 1.0, i = u,d, calculated 
by solution of Eq. l ]15| . 7 = 0.2 and k — in both plots. 
The cavity couplings Vi = qi — 1.0, i = A, B are chosen 
symmetrically. 

us to employ a quantum trajectory picture dividing the 
ensemble of quantum trajectories into two classes: The 
first class contains all realizations which have not been 
exposed to dephasing before the final photon emission 
whereas trajectories belonging to the second class have 
been dephased. The fingerprint of this unravelling in 
terms of quantum trajectories visible in the results for 
the ensemble average justifies the formally involved 
treatment of the open system's dynamics in terms of 
a PDP in the reduced system's Hilbert space. Fur- 
thermore, it has been shown that the long time decay 
of entanglement follows indeed an exponential decay 
which is a typical behaviour with Markovian dissipation. 
In this context the reciprocal proportionality of the 
measure of this decay to the dephasing strength n has 
been demonstrated. Asymmetries in both the coherent 
couplings of the electron-HH excitations to the cavity 
mode and the tunnel couplings between optical and 
lateral dots have been investigated without dephasing. 
The most important result of these simulations is that 
polarization-symmetry of the respective quantities is the 
crucial point concerning reliable production of entangled 
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photons. However, if this polarization-symmetry can 




Figure 17: V(t B ,t A ) (top) and P(t B ,t A ) (bottom) for s k u = 
1.1 / 3% = t\ = t k d = 1.0, k = A,B, calculated by solution 
of Eq. p§| . 7 = 0.2 and k = in both plots. The cavity 
couplings Vi = qi — 1.0, i = A, B are chosen symmetrically. 



not be assured, A-B-symmetry of the system's Hamil- 
tonian can also guarantee perfect entanglement if only 
simultaneously emitted photons are considered. 

In the future, it might be interesting as well as 
experimentally relevant to extend our analysis to the 
non- Markovian regime which would allow the modelling 
of further dissipation channels, e.g. hypcrfine interac- 
tion. The necessity of such a treatment depends on the 
particular dephasing channel that is the most relevant 
one of a given host material for the lateral and optical 
quantum dots. 



Acknowledgments 

We would like to thank Thomas Balder and Leo 
Kouwenhoven who contributed to early stages of this 
work and acknowledge financial support by the German 
Science Foundation (DFG). 



13 



1 M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

2 O. Benson, C. Santori, M. Pelton, and Y. Yamamoto, 
Phys. Rev. Lett. 84, 2513 (2000). 

3 O. Gywat, G. Burkard, and D. Loss, Phys. Rev. B 65, 
205329 (2002). 

4 T. M. Stace, G. J. Milburn, and C. H. W. Barnes, Phys. 
Rev. B 67, 085317 (2003). 

5 P. M. Visser, K. Allaart, and D. Lenstra, Phys. Rev. A 68, 
053805 (2003). 

6 F. Troiani, J. I. Perea, and C. Tejedor, Phys. Rev. B 74, 
235310 (2006). 

7 M. Larque, I. Robert-Philip, and A. Beveratos, Phys. Rev. 
A 77, 042118 (2008). 

8 F. Troiani and C. Tejedor, Phys. Rev. B 78, 155305 (2008). 

9 G. Pfanner, M. Seliger, and U. Hohenester, Phys. Rev. B 
78, 195410 (2008). 

10 K. Edamatsu, G. Oohata, R. Shimizu, and T. Itoh, Nature 
431, 167 (2004). 

11 S. M. Ulrich, M. Benyoucef, P. Michler, N. Baer, P. Gart- 
ner, F. Jahnke, M. Schwab, H. Kurtze, M. Bayer, S. Fafard, 
Z. Wasilewski, and A. Forchel, Phys. Rev. B 71, 235328 
(2005). 

12 R. M. Stevenson, R. J. Young, P. Atkinson, K. Cooper, 
D. A. Ritchie, and A. J. Shields, Nature 439, 179 (2006). 

13 R. J. Young, R. M. Stevenson, P. Atkinson, K. Cooper, 
D. A. Ritchie, and A. J. Shields, New Journal of Physics 
8, 29 (2006). 

14 A. Greilich, M. Schwab, T. Berstermann, T. Auer, R. Oul- 
ton, D. R. Yakovlev, M. Bayer, V. Stavarache, D. Reuter, 
and A. Wieck, Phys. Rev. B 73, 045323 (2006). 

15 N. Akopian, N. H. Lindner, E. Poem, Y. Berlatzky, J. 
Avron, D. Gershoni, B. D. Gerardot, and P. M. Petroff, 
Physical Review Letters 96, 130501 (2006). 

16 J. E. Avron, G. Bisker, D. Gershoni, N. H. Lindner, E. A. 
Meirom, and R. J. Warburton, Phys. Rev. Lett. 100, 
120501 (2008). 

17 R. J. Young, R. M. Stevenson, A. J. Hudson, C. A. Nicoll, 



D. A. Ritchie, and A. J. Shields, Phys. Rev. Lett. 102, 
030406 (2009). 

18 V. Cerletti, O. Gywat, and D. Loss, Phys. Rev. B 72, 
115316 (2005). 

19 M. Titov, B. Trauzettel, B. Michaelis, and C. W. J. 
Beenakker, New Journal of Physics 7, 186 (2005). 

20 C. Emary, B. Trauzettel, and C. W. J. Beenakker, Phys. 
Rev. Lett. 95, 127401 (2005). 

21 H. M. Wiseman, Ph.D. thesis, University of Queensland, 
1994. 

22 H.-P. Breuer and F. Petruccione, The theory of open quan- 
tum systems (Oxford University Press, Oxford, 2002). 

23 L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, Re- 
ports on Progress in Physics 64, 701 (2001). 

24 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Reviews of Modern Physics 
79, 1217 (2007). 

25 M. Blaauboer and D. P. DiVincenzo, Phys. Rev. Lett. 95, 
160402 (2005). 

26 C. Thelander, P. Agarwal, S. Brongersma, J. Eymery, L. 
Feiner, A. Forchel, M. Scheffler, W. Riess, B. Ohlsson, U. 
Gsele, and L. Samuelson, Materials Today 9, 28 (2006). 

27 R. Fiederling, M. Keim, G. Reuscher, W. Ossau, G. 
Schmidt, A. Waag, and L. W. Molenkamp, Nature 402, 
787 (1999). 

28 Y. Ohno, D. K. Young, B. Beschoten, F. Matsukura, H. 
Ohno, and D. D. Awschalom, Nature 402, 790 (1999). 

29 K. Blum, Density Matrix Theory and Applications 
(Springer, Berlin, 1996). 

30 C. W. Gardiner and P. Zoller, Quantum Noise: A 
Handbook of Markovian and Non-Markovian Quantum 
Stochastic Methods with Applications to Quantum Optics 
(Springer, Berlin Heidelberg, 2004). 

31 J. F. Clauser, M. A. Home, A. Shimony, and R. A. Holt, 
Phys. Rev. Lett. 23, 880 (1969). 

32 E. Jaynes and F. Cummings, Proceedings of the IEEE 51, 
89 (1963). 



